Long Range Interactions in Quantum Many Body Problem in One-Dimension: Ground 

State 
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We investigate the ground state properties of a family of iV-body systems in 1-dimension, trapped 
in a polynomial potential and having long range 2-body interaction in addition to the inverse square 
potential studied in the Calogero-Sutherland model (CSM). We show that for such a Hamiltonian, 
the ground state energy is similar to that of free fermions in a harmonic well with a displacement 
that depends on the number of particles and depth of the well. We obtain the ground state wave 
function and using random matrix results, study the particle density and pair correlation function 
(PCF). We observe that the particles are arranged in bands. Due to the presence of long range 
interaction, the PCF shows a departure from the CSM. 
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Theoretical understanding of the ground state proper- 
ties of complex many body systems have received consid- 
erable attention in recent years [10HSI!I1I1ES[I1 
ITU 0, El- In this context, we study rigorously a wide 
class of 1-dimensional TV-body systems having different 
densities, nature and strength of a 2-body interaction. 
We obtain the ground state properties of these systems 
having a 2-body potential V% = g /r\ 2 + $i , where m is 
the interparticle spacing and contains the long range 
2-body interaction. The system is trapped in a poly- 
nomial potential. This may be relevant in understand- 
ing the various aspects of the Bose-Einstein-Condensates, 
where a wide variety of potentials under a controlled en- 
vironment is possible. 

We derive the ground state eigenvalues and eigenfunc- 
tions for such systems and extract several interesting 
properties by identifying the square of the wave function 
with the joint probability distribution (j.p.d.) of eigen- 
values of non-Gaussian ensembles of random matrices. 
Using the polynomial method developed by Ghosh and 
Pandey in the context of random matrix theory (RMT) 
piU ITU , we observe band structure 0| in the particle 
density, which in turn corresponds to the density of ze- 
ros of the corresponding polynomials 0, 0, E| ■ For a 
given value of the interaction strength, we study the PCF 
for different interparticle spacings. Due to the presence 
of the longrange interaction, we observe a deviation from 
the CSM [1 0. 

We shall consider the ground state of a system of N 
particles satisfying the Schrodinger equation 



d 2 



N 



n 

i<j 



+ $i(xi,Xj) + Vi(xi) - E n 



where terms corresponding to long range interaction de- 



rive from the relation 



i¥=3 



PjXj) 
(Xi — Xj) 



(2) 



P{x) being any analytic function having a power expan- 
sion in x, while g and h are the interaction strengths. 
The system is trapped in a potential 



N 

V 1 (x l ) = Y / [P 2 (x i )-P'(x t )}. 

i=l 



(3) 



In this paper, we consider the case where P(x) is a poly- 
nomial of order 2m + 1 and is represented by 



(4) 



fc=0 
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For convenience, we take 02 m +i = 1- The parameter 7 
will determine the depth of the well. The case where 
m = corresponds to the CSM 1], where the 2-body 
interaction is purely of the inverse square type. For m > 
0, we will encounter the long range interaction. 

For such a Hamiltonian, the ground state wave func- 
tion can be written as 



-00 = 4>f = Y\_ \ x i ~ Xj\ X exp[- V" / P(U)dti}. (5) 

i<j i Jo 

It should be noted that if the particles are confined to 
a configuration space x\ > X2 ■ ■ ■ > a; at, the g/r 2 inter- 
^'actiB^rn one dimension does not allow them to cross, 
thereby respecting the ordering. This is valid provided 
the potential is not too attractive, which leads to the 
restriction g > -1/2 [TUTU. For such an ordering, the 
modulus becomes irrelevant and hence depending on the 
value of A (odd or even) the system is bosonic or fcrmionic 
in nature. Also, for such a configuration, ipo is nodeless 
(apart from the trivial ones at the points of coincidence) 



2 



and hence corresponds to the ground state of the system. 
Here 



= exp[- J2 / P{ti)dti\ 



(A 2 — A) = .g/2, 



and 



2A = -h. 
The ground state energy is given by 

E = ->yaiN(l + \(N-i)). 



(6) 



(7) 



(8) 



(9) 



(10) 



For ai > 0, this is the same as that of free fermions in a 
harmonic well 701, displaced by a factor r y\Nai(N — 1). 

First we prove that for P(xi) given by Eq.Q, Eq.JSJ) 
and Eq. ijTUI) are respectively the ground state eigenfunc- 
tion and eigenvalues. We study two cases, where we 
calculate tpo and £0 and recognize ip 2 as the j.p.d. of 
eigenvalues for the non-Gaussian ensembles of random 
matrices. Using the polynomial method, we study the 
particle density and the effect of long range interaction 
on the PCF. PCF shows a deviation from the CSM, which 
thereby suggests a departure from the "universal" result 
of RMT. 

The kinetic term of the Hamiltonian acting on 0>o gives 



-E 



dx 7 2 



> E 
- E 



dxi dx 
dx, 2 



(11) 



Substituting ipo from Eq. JSJ, the first term gives 



E 



d 2 {4>) 
dxi 2 



-[2(A 2 -A)$: 



i<3 



(Xi Xj ) 



The second term gives 



^E 



dxi ' dxi 



P{Xi) 



{xi Xj) 



while it is easy to see that the last term 



E 



dxi 2 



(13) 



(14) 



Thus we see that our choice of ipo diagonalises the Hamil- 
tonian (JTJ. Replacing P(xi) from Eq.Q in Eas. (jT3|) and 
(|14p . one can easily obtain the ground state energy for 
the cases where m > 0. 



We will consider the two cases, corresponding to 

P(x i )=~/(xf-a 1 Xi), 7>0, (15) 

and 

P(xi) = 7(a;f - a 3 xf + aiXi) j>0,a 3 <0. (16) 

They will not only illustrate the formation of multiple 
bands in the particle density for appropriate values of 
the parameter a\ , but also show the effect of long range 
interaction in determining the PCF. 

The Hamiltonian corresponding to Ea. (|15fl is, 



H 



2 



dxi 2 



n 

l<3 



<) 



(Xi 



E (^ ~ x rf 



i<j 



y 2 ^6 2fc2; f(17) 

k=l 



For 1 >> 7 > 0, the system behaves like an anisotropic 
oscillator, kept in a weak polynomial well. For such sys- 
tems, the motion is bounded for all physically possible 
energies. In the thermodynamic limit, i.e. for N — > 00 
and 7 — > 0, we will have a finite particle density, whose 
shape will depend on the well it is kept in. Finally, re- 
peating the steps outlined earlier, one can find that a 
choice of ipo 



■00 = 4>f = Y\ \ x * ~ x j\ x ex P 



1<J 



-7 E 




7 > 



(18) 

satisfies the Schrodinger equation, with g and h given by 
Eqs. JSJ and @ respectively, and 



b 6 = 1, 64 = -2ai, 62 = «i 2 - (6A + 3)/ 7 . 



(19) 



For 7ai 2 = 6A + 3 and 7 — ► 0, we can obtain the ground 
state wave function and eigenvalue for the central poten- 
tial V(ri2) = Ar\ 2 + B/r\ 2 and compare the result for 
N = 2 with that obtained in Ref . \Vq . For such a choice, 
we find Eq — 2(1 + X)y/ (37(2A + 1)) as compared to 
E = V^7(2A + 1) obtained in [l7|. The deviation is due 
to the effect of the potential well. For iV-particles, the 
ground state energy is given by Ea. i|l(J|) with a negative 
sign introduced due to the negative coefficient a\. The 
Hamiltonian corresponding to Eq. i|16fl can be written as, 



• : \ : 11 " 



dXl2 l <j( x i- x jf 



- a\ ( x i - x j) 4 



&2 ^2 (Xi - Xj) + «3 E ( Xi + X ^ + 72 E b 2k$$) 

i<j i<j 



k=l 



Here, the second, third and fourth terms have the same 
effect of making the particles remain in a bound state, 
with the parameter 7 and 03 controlling respectively the 
range and depth of the potential . 
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Repeating the steps outlined earlier, one can find that 
a choice of ipo 



V>o = JJkt - a;j| A exp 

i<j 



V — 
2-s 6 



EX; Xa 



with 7 > 0, 03 < satisfies the Schrodinger equation, 
with h given by Eq.J5J), and 



ai = A7/12,a 2 = — A 7 a 3 ,a3 = — 5A7/I2, 

610 = 1, b s = -2a 3 , 6 6 = (ag + 2a x ), 

fe 4 = -(20103 + (5 + 10A/3)/ 7 ), b 2 = a\+ a 3 (6A + 3^2/) 



The ground state energy is given by Ea.pOJI. As in the 
previous case, the parameter a\ will not only determine 
the position of the ground state but will also be crucial 
in determining the band structure. 

It is at this point that, we will write ipo in terms of the 
variables 



m=\- x 



l/(2m+2) 



X, j . 



(23) 



Then ipo can be written as 
^o = C 1 ^ 2 l[\y l ~y 1 \ x exp[-XY / Pi(U)dU}, (24) 



where 



P 1 (| /l )= 7 ^ C2fc+1 j / f +1 , 



(25) 



k=0 



is a monic polynomial of order 2m + 1, whose coefficients 
are related to those of P(xi) by 



C2/C+1 



(m-fc)/(m+l) 



0>2k+l- 



(26) 



Then ipo 2 is given by 

^o 2 = CjJlW - yjfexp[-/3j2 f % Pi{U)dti). (27) 



with /3 = 2A, and C being the normalization constant. 
Now one may interpret ■f/'o 2 to be identical with the j.p.d. 
of non-Gaussian ensembles of random matrices [15(. We 
define the n-particle correlation function 



Rn(yi,V2, 



- Un 



dy n+1 . . .dy N ip 2 . (28) 



as the probability of finding rt-particles in the intervals 
Hi and yi + Ay,, irrespective of the position of the other 
particles, n = 1 and 2 correspond to the particle density 
and PCF respectively. It is shown in Ref. 15] that R\(x) 



7 > Q(21) 
R 2 (r) 



- i 

■ n J 



f 

□ Ay =10"' 

• Ay=10" 3 

A Ay =10"' 

— GUE result 



FIG. 1: shows the smoothed PCF 01 f° r the Hamiltonian 

O, with P = 2 . N = 5 °. fl i = 60 and r = &V-Ri(y). The 
curves correspond to different values of Ay. The solid line 
corresponds to CSM (GUE) result, with Ay = 10" 2 . 



corresponds to the density of zeros of the polynomial hav- 
ing weight function given by 0. It has been shown by 
Dyson in the context of random matrices that for (3=1, 
2 and 4, R n can be written in terms of orthogonal and 
skew-orthogonal polynomials. For j3 = 2, the PCF can 
be written as 

N-l y 2 

-R2 (2/1,2/2) = V ^- 1 [g M (2/i)(Z M (2/2)]exp[-2 / P x (t)dt] 

(29) 

where (y) are orthogonal polynomials corresponding to 
the normalization condition 



q^(y)qv(y)Gxp[-2 Px{t)dt]dy = h^S^. (30) 
) Jo 

Now, we come back to the Hamiltonian discussed in 
Ea. (|17|) . For a\ > 0, the ground state energy is positive. 
For ai = a c = y / 2N/~/, we define E 0c = ja c N(l + \(N — 
1)). For Eq > Eq c , we observe the formation of two 
bands in particle density. Rescaling the result obtained 
by Pandey |16|], we get 

■kRM = 7MVW7 -(y 2 -ai) 2 - (31) 

For —Eoc < Eq < 0, a passage through the barrier 
is possible resulting in a splitting of each of these levels 
into two neighboring ones, corresponding to the state in 
which the particles move simultaneously through both 
the barriers. This corresponds to a single band case with 
a dip around the origin. Finally, for Eq < —Eq c , we get 
a single band with a maximum at y = 0. This is given 
by 



— Wai 2 + 6N/ 1 -2a 1 } + y 2 
3 7 



— {Vai 2 +6A7 7 + ai}-2/ 5 
3 7 



1/2 



,(32) 
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FIG. 2: Particle density for the Hamiltonian given in Eq.(20), 
with N — 50. For a given as and different values of the 
parameter a\ , we observe a shift from single band to multiple 
band structure. 



for ai < y / 2N/j for the single band case. We observe 
that for finite R\ (y) , letting N — > oo implies that 7^0 
as TV -1 . It should be noted that in the thermodynamic 
limit, the particle density, after proper scaling, is inde- 
pendent of (3 and hence applicable for both bosons and 
fermions. However the calculation of the PCF, essential 
for the study of thermodynamic properties, is restricted 
to (3 = 1, 2 and 4. It has been proved in Ref.0,0] that 
for interparticle spacing (yi — y%) = Ay — > 0, N — > 00 
and defining r = Ay.Ri(y), the scaled PCF for CSM is 
universal. For (3 = 2, it is given by 



Mr) 



sin 2 (7rr) 
(nr) 2 



(33) 



For these systems with Ay = 10 -7 , we are in the re- 
gion of the spectrum where the particle density is high 
and hence interparticle spacing low. Here, PCF obeys 
Ea. (|33|l . However, for Ay = 10~ 3 and 10 _1 , we are in 
the tail region of the two bands. It is here that we ob- 
serve a distinct deviation from that observed in CSM. 



This is due to the contribution of the long range interac- 
tion, where the third term in Ea. (|17[) contributes to the 
PCF. This becomes negligible as one makes Ay smaller, 
until it disappears for Ay = 10~ 7 . 

For the Hamiltonian corresponding Ea. (|20() . it can be 
shown analytically that for large N, that the particle 
density can be written in terms of the higher moments 
as 



%Ri(x) = 7(-x 10 + 2a 3 a; 8 - (a 3 2 + 2ai)x 6 

+ (2iV/ 7 + 2a 1 a 3 )x i - {a x 2 - 2M 2 / 7 + 2a i N/ 1 )x 2 



+ {2 ai N + 2M± - 2a 3 Af 2 )/ 7 ) 1 / 2 

where M p = J x p Ri(x)dx, x being the scaled variable. 
These moments can be calculated numerically to obtain 
Ri(x). However, a more convenient way is to directly 
construct the polynomials through the recursion relation 
and to do the sum of Eq.J2H|), with x — y. For a 3 < 0, we 
define a\ = a c = a 3 /4, which corresponds to a ground 
state energy E 0c = -ja c N(l + X(N - 1)). For E < 
Eq c , we observe a single band structure in the particle 
density. For Eq > Eo c , (i.e. a% < a c ), the particles 
experience a repulsion at the center and thus we see a 
transition from single band to three band structure. As 
Eq << Eq c , particles are completely repelled from the 
center and finally collect at the two end of the spectrum, 
thereby giving rise to two bands. 

Thus we have studied, for both bosons and fermions, 
the ground state properties of a class of iV-body systems 
with long range two body interaction in addition to the 
inverse square potential. We observe a distinct deviation 
from the "universal" result for the PCF, as obtained in 
the CSM. We also observe that particles are arranged 
in bands. To study different thermodynamic properties, 
it is necessary to study the excitation spectrum of such 
Hamiltonian. At hnite temperature and finite N, we ex- 
pect similar smoothing of the oscillations in the particle 
density 0] as observed in Ref.0|- We will come back 
to this in a later publication. 

I am grateful to Dr. P. K. Panigrahi for many useful 
discussions. 
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